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ABSTRACT 

In  a  recent  paper*  by  T.  J.  Harvey,  there  is  presented  the  idea  of 
employing  known  solutions  of  free  vibrations  of  nonlinear  single-degree 
conservative  systems  to  obtain  solutions  of  forced  vibrations  of  these 
systems.   The  forcing  function  is  taken  to  be  of  the  same  form  as  the 
free  response  and  the  solutions  are  obtained  quite  easily.   Harvey  makes 
the  point  that  while  this  excitation  is  not  the  sinusoidal  excitation 
that  many  analysts  have  attempted  to  study,  for  many  purposes  there  is 
reason  to  prefer  some  alternate  form  of  excitation  if  by  doing  so  useful 
solutions  are  readily  obtained.   In  a  paper  of  several  years  back,  M.  A. 
Sadowsky*  derived  the  characteristics  of  several  nonlinear  spring  systems 
for  which  the  free  vibrations  could  be  obtained  in  rather  simple  closed 
forms.   In  particular,  Sadowsky  introduced  two  new  springs  which  give  solu- 
tions in  terms  of  the  elementary  functions  and  several  others  which  give 
solutions  in  terms  of  the  Jacobi  elliptic  functions. 

Since  there  are  actually  only  a  few  fully  developed  cases  of  forced 
nonlinear  vibrations  available  to  analysts  for  testing  the  power  and  ac- 
curacy of  analytical  methods  and  for  other  purposes,  it  appeared  to  be  of 
interest  to  apply  the  viewpoint  of  Harvey  to  several  of  the  particular 
cases  given  by  Sadowsky  so  as  to  work  out  all  the  details  of  the  corres- 
ponding cases  of  forced  vibration  and  thus  afford  several  new  examples. 
The  major  content  of  this  thesis  is  the  full  development  of  these  cases. 

In  addition,  however,  the  viewpoint  of  Harvey  is  applied  to  two  cases 
of  piecewise  linear  springs  simply  for  the  purpose  of  illustrating  such 


*See  references  on  page  41. 
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applications  and  for  the  purpose,  as  well,  of  obtaining  additional  cases 
of  nonlinear  vibrations. 

Also,  since  the  forcing  functions  differ  from  the  sinusoidal  func- 
tions that  many  analysts  have  used,  it  was  considered  desirable  to  obtain 
and  exhibit  a  simple  measure  of  the  degree  of  "distortion"  of  the  actual 
forcing  function  as  compared  to  the  sinusoidal  shape.   Likewise,  the  dis- 
placements are  not  sinusoidal  functions  of  time,  and  the  "distortion"  of 
the  response  (as  compared  to  a  sine  curve)  is  also  examined. 

This  work  was  done  largely  during  the  spring  of  1961  at  the  U.  S. 
Naval  Postgraduate  School,  Monterey,  California.  The  writer  wishes  to 
express  his  appreciation  to  Professor  J.  E.  Brock  for  having  suggested 
the  topic  and  for  guidance  during  its  development. 
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SYMBOLS 

t  *  time  a  independent  variable,  chosen  so  that  $(<?)  =  0 

0  ,~c  s  diraensionless  independent  variables 

i^r)  s  deflection,  normalized 

p(t)  ■  external  force  applied  to  the  system 

5 1})    =  spring  function  such  that  S(o)-   0  and  $(o)  =  1.0 
$(#■)    has  the  same  dimension  as  does  :? 

f>     -   assignable  dimensionless  constant  as  defined  in  the  relation  (5) 

ttl    =  mass  or  moment  of  inertia  of  the  system 

4{     =  spring  constant  for  a  small  deflection 

(lit    ss  l^/m    =  natural  frequency  for  a  small  free  oscillation 

7^  ■  ~*"J''   ~     period  of  a  small  free  oscillation 

To  =  period  of  a  free  vibration  (  P  =  0) 

U/    a  frequency  of  vibrations 

7*  *  -Sil-  s  period  of  forced  vibrations 

J  «yflGjj  ■  amplitude  of  deflection 

—  denotes  the  evaluation  at  "jt -   % 

o(  =  modulus  of  Jacobi  elliptic  func*  ion  when  a  solution  contains 
such  functions  (refer  to  Table  1) 

o(,  ss  modulus  of  Jacobi  elliptic  functions  whe^  a  spring  function 
contains  such  functions  (refer  to  Table  1) 

K    -   the  complete  elliptic  integral  of  the  first  kind  with  modulus  <V 

Kj  -   the  complete  elliptic  integral  of  the  first  kind  with  modulus  o<; 

/|  ■  the  distortion  of  fit)   wave  form  as  defined  in  Eq.  (17) 

B  ss  the  distortion  of  f$t)    wave  form  as  defined  in  Eq.  (18) 

8    as  a  constant  for  the  springs  of  cubic  parabola  types  such  as 

RYj  -Y  +  PY* 

Y  ~   deflection  which  is  not  normalized 
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1.   Introduction. 

*  - 

The  equation    7t1  f    -f  k  $    -  Fc  C^ft^t  (1) 

describes  the  behavior  of  a  linear  single-degree-of- freedom  conservative 
system  subjected  to  a  sinusoidal  excitation. 

A  natural  extension  for  the  study  of  a  nonlinear  conservative  system 
is 

m  jf  +^Sft)  =  Ji  <u*k)t  (2) 

where      is  a  nonlinear  spring  function,  such  as   S  .   Many 

investigators  have  studied  Eq.  (2)  in  detail  and  obtained  solutions  which 
indicate  the  effects  of  parameter  U)   upon  %■     .   But  Eq.  (2)  is  complicated 
to  solve  analytically,  and  in  most  cases  the  solutions  are  obtained  by 
approximation  methods. 

T.  J.  Harvey  (ref.  1)  suggested  a  simple  analytical  method  of  study- 
ing the  forced  vibrations  of  nonlinear  systems  by  making  use  of  the  solu- 
tions for  free  vibrations,  i.e.,  where  the  right  hand  side  of  Eq.  (2)  is 
zero.   Some  years  ago  M.  A.  Sadowsky  (ref.  2)  presented  simple  solutions 
of  free  vibrations  for  various  nonlinear  springs. 

This  thesis  applies  Harvey's  method  to  several  of  the  cases  given  by 
Sadowsky,  thus  obtaining  several  interesting  solutions  for  forced  vibra- 
tion of  nonlinear  systems. 

In  addition,  a  study  is  done  on  the  following: 

(a)  The  wave  form  of  $(t)    is  compared  with  sinusoidal  one  which 
is  that  of  a  linear  system. 

(b)  Traditionally  the  right  hand  side  of  Eq.  (2)  is  a  sinusoidal 
function,  as  shown.   By  Harvey's  method  it  differs  from  a  sinusoidal  func- 
tion, but  is  still  a  periodic  function.   These  differences  are  compared. 
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(c)   Harvey's  method  can  be  applied  to  other  than  Sadowsky's 
nonlinear  springs.   It  is  applied  to  the  spring  characteristics  which 
consist  of  linear  segments,  as  examples. 


2.   B^sic  equations. 

The  forced  vibration  of  an  undamped  nonlinear  single-degree-of- free- 
dom system  is  described  by  the  equation 

'i  +  ks(t)  =  fit)  o) 

where  5f#.J  *-s  a  nonlinear  spring  function  normalized  so  that  $(0)  ~  0 
and  $\o) «•  \*d    »  a^d  Pfx)  is  the  external  force  applied  to  the  system. 
Eq.  (3)  is  simplified  by  substituting: 

^  ^   Mi    and    T.  (4) 

Then  JfL    +  ,S(7)  ^  XFW  (3a) 

Eq.  (3a)  as  it  stands,  is  of  great  difficulty.   Even  if  p(r)  is 
periodic  and  if  we  presume  (by  way  of  analogvwith  linear  systems)  that 
slight  damping  present  in  the  system  but  not  accounted  for  in  the  mathe- 
matical description  causes  the  "transient'1  to  die  outs  the  "steady  state* 
solution  is  generally  not  obtainable  except  by  difficult  analyses  leading 
to  approximations.   However,  the  analysis  is  greatly  simplified  by  henceforth 
confining  attention  to  only  such  f(f)     as  satisfy  the  relation 

i  *  "  A  sew]  (5) 

where  p     is  a  dimensionless  constant,  either  positive  or  negative.   Such 
a  forcing  function  is  called  a  "natural  forcing  function"  (ref.  1).   A 
positive  p  means  the  external  force  is  in-phase  with  restoring  force,  a 
negative  P     means  out-of-phase  (180),  and  ^-  0    means  free  vibration. 
Substituting  the  relation  (5)  into  Eq.  (3a),  we  obtain 

"JTT  +0-f)Wj=  O  (6) 

This  is  further  simplified  by  putting 

G    «  //-(°  "C  (7) 


.ft 

Then        -jfx     f  3#)  -  (8) 

which  is  described  as  the  case  of  the  pseudo-free  oscillation.   To  solve 
Eq.  (8),  let  the  initial  conditions  be 

The  potential  energy  stored  in  the  spring  is 

Uft)=  \Sty)dfr    t        U(o)  =C 

•So 
Furthermore ,  by  conservation  of  energy,  Vc      is  related  as 

v;a  =  z  u(j) 

where  .'?  is  the  amplitude  of  deflection.   If  T/4  is  defined  as  the  quarter- 
period  of  oscillation,  H-    can  be  written  as   ft    -  ftCJfa.) 
Eq.  (8)  is  integrated  and  simplified  as 

Integration  of  Eq.  (9)  leads  to 

ft     df 

&        -=  I    [2  «?)-*«(*)}*  <10) 

Again  Eq.  (10)  can  be  written  in  the  following  general  form  with  ]Z 

as  a  parameter, 

&        —  9($,    J)  (Uy 

The  inverse  form  of  Eq.  (11)  can  be  expressed  as 

Since  the  dimensionless  variable  Q     is  closely  related  to  time  ts 
using  Eqs.  (4)  and  (7), 

e  =u)*t\i-f>  da) 

Then  #  =  tH+Ff*  i)  (13) 


This  is  the  solution  of  the  forced  vibration  of  Eq.  (3)  in  which 
external  forces  f(t)    satisfy  the  relation  (5).   For  free  vibrations  fit)-': 
and       .   Then  the  solutions  of  free  vibrations  are  from  Eq.  (13), 

*-*(**..?;  (U) 
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Comparing  Eqs,  (13)  and  (14),  the  solutions  of  forced  vibrations  are 
obtained  from  those  of  free  vibrations  by  multiplying  the  independent  vari- 


able t(^  time)  by  the  factor  \fT~~P      •   The  assignable  dimensionless  con- 
stant t^      is  the  ratio  of  the  external  force  to  the  restoring  force,  namely 

where  the  relation  (5)  is  rewritten.   The  amplitude  of  external  force  p 
is  given  as  2.0^  or  3.5^  ,  for  example,  which  means  the  force  unit  corres- 
ponding to  the  value  of  the  spring  constant  (at  small  deflection)  multiplied 
by  2.0  or  3.5  units  of  deflection  respectively.   Or  in  other  terms,  Eq.  (3a) 
represents  a  forced  vibration  wherein  S(f)  plays  the  role  of  a  restoring 
"force"  and  ~%  \-(X)  plays  the  role  of  an  exciting  force;  the  time  scale 

has  been  chosen  so  that  both  "forces"  have  the  same  dimension  as  the  "de- 
flection response",  frfi).     Thus  we  can  say  that  the  amplitude  of  exciting 
force  (  F/%^       )  has  magnitude  2.0  or  3.5,  say.   For  every  amplitude  of  de- 
flection ^   ,  the  values  of  o— •$(?/   are  found,  which,  in  turn,  deter- 
mines the  assignable  constant  P    . 

Sadowsky's  paper  gives  the  essential  details  for  ten  cases  of  nonlinear 
springs  for  which  the  free  oscillations  are  relatively  simple.   Two  of  the 
solutions  are  in  terms  of  elementary  functions,  and  the  others  involve  the 
Jacobi  elliptic  functions.   Corresponding  to  each  of  these  cases,  by  using 
Harvey's  method,  it  is  possible  to  construct  a  complete  solution  for  the 
forced  oscillations  resulting  from  "natural"  excitation.   In  this  thesis 
four  of  the  cases  are  developed  in  detail--the  two  cases  involving  only 
elementary  functions  and  two  of  the  cases  (cubically  hardening  and  cubi- 
cally  softening  --  i.e.,  Duffing' s  spring)  involving  elliptic  functions. 
Also,  two  cases  of  piecewise-linear  springs  (cases  not  included  among  those 
given  by  Sadowsky)  are  also  developed. 
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Sadowsky's  solutions  are  expressed  in  terms  of  the  initial  velocity 
\/c        ,  which  can  be  replaced  by  Y       ,  for  convenience,  with  the  energy  reU 
tion  V0*  —  -2  ^(?j   in  the  conservative  systems.   The  solutions  are  tabu- 
lated in  Table  1  and  their  approximate  spring  characteristics  are  shown 
in  Fig,  lo   The  symbol  o(      represents  the  modulus  of  Jacobi  functions, ft 
means  the  complete  elliptic  integral  of  the  first  kind  with  modulus  o< 

In  Fig.  1,  the  numbers  in  circles  represent  the  spring  numbers  as  in- 
dicated in  Table  .   Springs  marked  as  1,  2,  3,  and  4  are  drawn  heavily. 
The  family  springs  7  is  in  between  springs  1  and  4  depending  upon  the  modu- 
lus o(i      »  where  for  o(\-0    it  coincides  with  spring  1  and  for  o(,=  |,0 
it  coincides  with  spring  4.   The  space  between  springs  4  and  3  is  filled 
with  the  family  of  springs  6  depending  upon  the  modulus  ofi   .   For  o'i-O 
it  coincides  with  spring  3,  and  for  c(r-l'0    it  coincides  with  spring  4. 
Each  one  reaches  a  maximum  and  drops  down  to  zero.  The  family  of  springs 
5  is  of  softening  type,  and  so  is  the  elementary  spring  2  obtainable  with 
>Ai  ~l<0      .   Springs  8,  9,  and  10  are  cubic  parabolas,  they  are  not  shown 
in  Fig.  1  since  they  can  readily  be  visualized. 


Table  1  -  The  solutions  of  forced  vibrations  for  Sadowsky's  springs 
(adapted  from  Sadowsky;  ref.  2) 

The  solutions  of  forced  vibrations  such  as 

m  j   f  k  Stt)  -   Fb) 

are  tabulated  below,  in  which  F(t)      satisfies  the  relation  (5). 


S(Z) 
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Fig.  1  Characteristics  of  Sadowsky's  springs, 
(from  Sadowsky,  ref.  2). 


i.o- ■ 


to- 


Fig.  2  Distortion  of  the 
wave  form.  A  is 
defined  in  Eq.  (17). 


Fig.  3  Distortion  of  the 
wave  form.   B  is 
defined  in  Eq.  (18). 


A.   Frequency  response. 

When  the  exciting  force  is  such  that  the  relation  (5)  is  satisfied, 
the  period  of  forced  vibration  is  readily  obtained  from  the  solution  of 
free  vibration.   If  the  period  of  free  vibration  is  obtained  as: 

T.  =  Ut) 


then  by  comparing  Eqs.  (13)  and  (14),  we  see  that 

T       I 


(15) 


Since  To  is  varying  with   ^  ,  the  ratio  {^fj    may  be  related  to  the 


period  of  a  small  free  oscillation  Tn 

Tr,         \T*/   j~ 


n 

i 


(16) 


This  means  that  corresponding  to  a  single  frequency  response  curve 
of  if   versus  (— r-J   ,  obtained  for  free  vibrations,  we  obtain  a  full 
family  of  response  curves  for  forced  vibrations,  i.e.,  curves  of  ^versus 
[ X.  J   for  various  values  of  P 


B.   Wave  form  of  f(t)    . 

When  the  exciting  force  is  such  that  the  relation  (5)  is  satisfied^ 
the  solution  is  exactly  the  same  as  that  for  free  vibration  except  for 
the  change  in  frequency,  which  is  >J|- f>        times  that  of  the  free  vibra- 
tion.  Therefore  the  wave  form  of  jf^r)  depends  upon  the  amplitude  of 
deflection  ^  only,  and  not  upon  the  amplitude  of  exciting  force.   When 
if    is  small,  the  springs  are  very  nearly  linear,  and  the  wave  forms 
are  almost  purely  sinusoidal.   As  ^    increases,  the  wave  forms  of  %(t) 
are  distorted  from  the  sinusoidal  forms. 

Furthermore  when  the  spring  characteristics  are  symmetrical,  the 
wave  forms  are  quarter-wave  symmetrical.   In  addition,  they  tend  to  lie  be- 
tween a  rectangular  and  a  triangular  wave  form. 

To  show  the  distortion  of  <j{t)  wave   forms,    a  measure  A   is   defined  as 

r\  ~  ^[tJJc    -jt  dT  (i?) 

If  we  let  the  quarter  period  T/4  correspond  to   ' "   ,  the  measure  A 


represents  an  area  under  the  curve  of  -~A     as  shown  in  Fig.  2.   Note 
that  A— /  »  f°r  a  linear  spring  or  a  small  oscillation. 
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C.   Wave  forms  of  f- (i~)    . 

When  the  external  force  f-(t)    satisfied  the  relation  (5),  the  wave 
form  of  \-Lt)     is  the  same  as  that  of  the  spring  force  51^1    .   Therefore 
it  depends  on  #   only,  not  the  amplitude  of  the  external  force  p 
For  a  symmetrical  spring,  the  wave  form  of  PCTj  is  quarter-wave  symmetri- 
cal in  the  same  manner  as  that  of  ^  It)   .   For  a  small  oscillation  fft) 
is  almost  a  pure  sinusoidal  form,  from  which  the  wave  form  is  distorted 
as  ft       increases. 

To  compare  this  wave  form  with  the  traditional  forcing  function  of 
sinusoidal  form,  a  measure  B  is  defined  as 

b  -imp -f1^ 

<4    F  (18) 

Note  that :  Q—l-D    for  a  linear  spring  or  a  small  oscillation. 
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3.   Examples, 

The  results  of  six  examples  are  presented  in  Fig.  (4)  through  Fig,  (15) 
Figures  4,  6,  8,  10,  12,  and  14  are  all  similar  representations  as  follows; 

(a)  Ordinates  are /f  (or  ^/c      where  C  is  a  constant). 

(b)  Solid  curves  represent  frequency  response  curves  for 
labled  values  of  p «  FO^t)     >  including  f  =■  0      »  the  case 
of  free  vibration. 

In-phase  forced  vibrations  are  indicated  by  positive  values  of  p~      , 
out-of-phase  by  negative  values.   These  curves  relate  %      (ordinate)  to 
s^i\       (bottom  scale  of  abscissa). 

(c)  Distortion  of  $(t)  is  indicated  by  a  measure  A,  distortion 
of  p(t)  by  a  measure  B,  both  functions  of  ^   alone,  and 
plotted  on  vertical  scales. 

(d)  The  dotted  curves  show  spring  characteristics  5(<f)  (top 
scale  of  abscissas)  as  functions  of  ^   .   In  Figs.  12 
and  14,  this  relation  is  shown  to  smaller  scale  in  the 
inset  and  not  in  full  scale  as  in  the  preceding  figures. 

Figures  5,  7,  9,  11,  13,  and  15  actually  display  the  distorted  wave 
forms  for  ft)    and  f(t)   . 
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A,   Example  1     ${f)  =  ^J 

The  characteristic  of  the  spring  is  extremely  hardening  as  shown  in 
Fig.  4  (dotted  curve).   The  solution  of  the  forced  vibration  is  already 
tabulated  in  Table  1,  namely 

jmf^   A^i}  A^  (19) 

where 

W   =  (4i  A&t  }')  JT^"  ,    ^  =  JL.  (20) 

(a)  Frequency  response. 

From  Eq. (:0)  we  can  get  the  solution  which  indicates  the 
effects  of  frequency  upon  !f-   ,  making  use  of  T*  ~uT      and  T*  -  ~^f    .   In 

the  dimensionless  form,  we  get 

/ 
—  (21) 

For  the  given  value  of  p   ,  P    is  obtained  by  dividing 
p   by  S(H)    »  then  the  ratio  J£   is  obtained  from  Eq.  (21).   The  fre- 
quency response  curves  are  plotted  in  Fig.  4. 

(b)  Wave  form  of  !f(t)     . 


(22) 


Eq.  (19)  is  rewritten  in  a  convenient  form, 
%(i)  aj*C'  (  ***  Ja*h  tut) 

For  a  small  oscillation,   M*«K    -4*P-  »  *+n  tit 
With  Eq.  (22)  the  wave  forms  of  ,f(t)      are  plotted  for  §  ■  0, 
0.5,  1.0,  1.30,  and  J£  in  Fig.  5(a).   Since  Eq..  (22)  shows  that  %tt)      is 
quarter-wave  symmetrical,  plotting  is  done  for  a  quarter  period.   For  a 
small  oscillation,  it  is  clearly  a  pure  sinusoidal  form,  and  as        in- 
creases it  is  distorted.   To  show  the  distortion,  the  measure,  A(as  defined 
in  Eq.  17)  is  found.   The  integral  is  not  easily  performed  analytically 
so  that  a  numerical  integration  method  was  used.   The  measure  A  is  a 
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function  of  ^f-  ,  and  is  plotted  along  ^  in  Fig.  4. 

(c)  Wave  form  of  f;t) 

Using  Eq.  (19)  and  the  relation  (5), 

rearranging  these  into  a  dimensionless  form,  we  get 

fit)  a**-  «!t 

~F~~  Q+ to} equity*  (23) 

For  a  small  oscillation,  ^%a     &  '    '   =  ***  <dt 

q  < 

which  is  purely  sinusoidal.   As  4.  increases,  it  is  distorted.   With  Eq. 
(23)  the  wave  form  of  f(t)    is  plotted  in  Fig.  5(b). 

To  compare  fit)   wave  form  with  the  traditional  sinusoidal  function, 
the  measure  B(as  defined  in  Eq.  18)  is  obtained: 

which  is  a  function  of  %      only,  and  is  plotted  in  Fig.  4. 

(d)  Discussion. 

This  spring  has  an  extremely  hardening  characteristic.  At 
ft  —  •£-   ,  the  spring  force  becomes  infinity,  its  slope  is  infinity,  and 
the  measure  B  is  zero.   For  slightly  smaller  values  of  H.        ,  the  me*  sure 

Q 

B  is  far  away  from  unity.   Fig.  5(b)  shows  that  the  wave  form  of  p(t)  may 
become  quite  distorted. 

When  ft    is  far  away  from  the  value  of  "tyz    ,  the  measure  B  is  not  far 
away  from  unity,  and  Fig.  5(b)  shows  the  wave  form  of  fit)   is  not  so  much 
distorted. 
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B.   Example  2„     6^J  =  "^TT" 

This  spring  has  a  softening  characteristic  as  shown  in  Fig.  6, 
The  solution  of  the  forced  vibration  is  already  tabulated  in  Table  1. 

■  *f   ^  -link  J  Ami  u/t  ^24) 

where 

u/^(<4.44?)/7T7    >     p=?  (25) 

&  S 

(a)  Frequency  response. 

From  Eq.  (25)     \   =  Co4fj=L==r  (26) 

For  a  given  value  of  p  ,  P  is  obtained  as  a  function  of 
"Q  ,  and  the  dimensionless  values  of  pr  are  found.  The  frequency  re- 
sponse curves  are  plotted  in  Fig.  6  forp=tf  and  i&|^. 

(b)  Wave  for  of  }f(t)  . 

Eq.  (24)  is  rewritten  in  a  convenient  form, 

For  small  oscillation,   "*g!^  -#-^  =  4**  *]T 
which  is  purely  sinusoidal.   With  Eq.  (27)  the  wave  form  of  if(t)   is 
plotted  in  Fig.  7(a)  for  f**&   ,  1.0,  2.0,  5.0,  and  00,  As  ft      increases 
it  becomes  distorted.  And  as  ft-    approaches  infinity  it  goes  to  a  rec- 
tangular form. 

To  show  the  distortion,  the  measure  A  is  found  and  plotted  along  b- 
in  Fig.  6.   The  integration  to  find  the  measure  A  is  performed  by  a 
numerical  method. 

(c)  Wave  form  of  f(f)   . 
Using  Eq.  (24)  and  the  relation  (5),  we  get 

F(t)m  Pks{SW]  ,     p.  pksm 


rearranging  these  into  a  dimensionless  form, 
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(28) 


For  a  small  oscillation,    *r*  .  -=~  -  ^  *#" 
which  is  purely  sinusoidal.   As  Ff    increases  it  is  distorted.   With  Eq. 
(28)  the  wave  form  is  plotted  for  f   =  0,  0.5,  0.66,  and  1.0  in  Fig.  7(b). 

To  compare  f(f)   wave  form  with  a  sinusoidal  one,  the  measure  B 
is  found, 

3  =  u4f  (29) 

which  is  plotted  along  ^    in  Fig.  6, 

(d)  Discussion. 

The  spring  is  of  softening  type;   at  %■    »  0.66  5$)  reaches 
its  maximum  value  (which  is  equal  to  0,384).   Beyond  this  point  it  decreases 
and  the  spring  force  approaches  zero  as  £  approaches  infinity  as  shown  in 
Fig.  6.   Up  to  tf     =  0.66,  the  measure  B  is  not  far  away  from  unity,  and 
Fig.  7(b)  shows  that  the  wave  form  of  p(t)   is  not  distorted  appreciably. 

Far  beyond  the  point  ^  =  0.66,  the  value  of  B  is  much  greater  than 
unity  and  the  wave  form  is  greatly  distorted. 
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Co   Example  3.   Duffingcs  equation.    S($)   ^  (f   +  "^  $- 

When  the  spring  characteristic  is  a  cubic  parabola,  the  equation  of 
forced  vibration  is 

™r  +  ktf +(!(*)  ■=  FW  <30> 

in  which  Eq.  (3)  is  rewritten  here.  X    *s  tne  deflection.   If  it  is  ex- 
pressed in  inches 9  6  is  in  inch  ,  and  -k     in  lbs/ in. 

To  solve  Eq.  (30)  with  the  right-hand  side  set  equal  to  zeros  the 
general  cubic  parabola  (Y'+PT  J  can  De  expressed  in  the  form  of  Table  1, 
i.e. ,  (  #  +  3-  %  )      by  proper  substitution  regardless  of  ft    . 

ut       Y=7tt7*  (31) 

Then  the  equation  of  free  vibration  becomes 

in  which  |  is  a  dimensionless  quantity. 

Then  we  obtain  the  solution  of  forced  vibration  immediately  from 
Table  1. 

p)  =  }c«(tot-K,ol)  C33) 

where       (V   «  4  (  '  +x  ^TV*3^" 

o<  =.  £(lt±-jj-  y       as  the  modulus  of  Jacobi  functions      (34) 
f  —  ~^=r  {Tp  *-n  a  dimensionless  form 

And  K  is  the  complete  elliptic  integral  of  the  first  kind  with  modulus 
o(   ,  i.e.,    «  =    (i-ofV*,^)4  j$ 
(a)  Frequency  response 

From  Eqs.  (33)  and  (34)  with  J"=  -^-  ,  we  get 


In  Fig.  8  the  frequency  response  curves  are  plotted  for  f  =  c    »  A  ^feT? 
and  A  2  ^  g 
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(35) 


Wave  form  of  %(£)  • 
Eq.  (33)  is  rewritten; 

For  a  small  oscillation,     .c  ~j        '-  ^"^ 
which  is  purely  sinusoidal.   As  $■      increases  it  becomes  distorted. 

However  there  is  a  limit  in  the  distortion.  As ^  approaches  infinity  the 

i 
value  of  o<   does  not  exceed  f£      .   In  Fig,  9(a)  the  wave  form  of  ffi) 

is  plotted  for  #  ■  0  and   fr    =  Go  . 

To  show  the  distortion  the  measure  A  is  founds  /|  =  fjL) — ~r-— 
which  is  plotted  in  Fig.  8  along  H-    .   It  does  not  change  appreciably 
from  unity.   Even  if  ^  approaches  infinity,  A  approaches  0o942„   In  other 
words,  the  wave  form  of  f(t)     is  not  distorted  much  from  the  pure  sinusoidal 
form. 

(c)  Wave  form  of  ffa) . 

From  Eq.  (33)  and  the  relation  (5), 

rearranging  these  into  a  dimensionless  form,  we  get 

-£&1  «    '     f^^-Kj  +  iF^^-Wl       (36) 
with  the  modulus  q(     „ 

For  a  small  oscillation,  ^t%d  p  A^fi^T  which  is  purely  sinu- 
soidal. With  Eq.  (36)  the  wave  form  of  f(t)  is  plotted  for  J  -  0S  1.0, 
2.0,  and  infinity  in  Fig.  9(b). 

To  compare  f(t)   with  a  sinusoidal  force,  the  measure  B  is  found, 


r   _  -X  (Jl±±X- 


which  is  plotted  along  ft     in  Fig.  8. 
(d)  Discussion. 

This  spring  is  a  hardening  type  as  shown  in  Fig.  8  (dotted 
curve).   Every  type  of  cubic  parabola,  such  as  J+ftf   ,  can  be  reduced  to 
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3  J 

the    form    5/^j  =  f  +  2  ft  by   letting     \   ~{=nf  fr      without    losing  the 

generality  of  the  governing  equation    (Eq.    3). 

Fig.  9(b)  shows  that  the  wave  form  of  f(t)    is  not  distorted 
to  an  appreciable  amount,  and  the  measure  B  is  not  far  away  from  unity. 

The  application  of  Harvey's  method  to  Duffing' s  equation 
was  considered  in  a  recent  note  (ref.  3)  which  appeared  after  this  thesis 
was  nearly  completed.   For  the  hardening  spring,  considered  in  this  sections 
the  results  given  in  the  note  agree  with  those  given  here.   However,  in  the 
case  of  the  softening  spring,  treated  in  the  next  section  following,  the  re- 
sults do  not  agree  since  the  note  contains  an  error.   In  place  of  the  en 
function  in  the  note,  the  cd  function  should  appear  for  the  softening  case8 
and  the  argument  of  the  function  should  also  be  slightly  modified. 
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Dc      Example   4,      Duffing's   equation.        S[$)—   $  ~  2  ft 

When  the  spring  characteristic  is  a  softening  type  of  cubic  parabola, 
the  general  characteristic  equation  is  expressed  as 

Rm  =  Y-  p  r1 

which  can  be  reduced  the  form  of  Table  1  by  letting  J  =  f^-fi  ft  ,  where 
<f  is  a  dimensionless  quantity.  Then  without  losing  the  generality  of 
the  governing  equation  (Eq.  3),  we  get 

The  solution  of  forced  vibration  is  already  tabulated  in  Table  1, 

%(t)    =  }  M(uJtM)  (37) 


where 


rr^X, 


*=}(*■-?)'*  (38) 

(a)  Frequency  response. 

From  Eq.  (37)  and  (38)  with  7=  -^-,  we  get 

-^- mirror)  (rr=) 

where  K  is  the  complete  elliptic  integral  of  the  first  kind„   The  fre- 
quency response  curves  are  shown  in  Fig.  10  for  p  -  0  and  ±  0.2  ytfzp 

(b)  Wave  form  of  ffi)   • 

Eq.    (37)    is   rewritten   in  a  convenient    form^ 
-W~=   M(WtM)  where     o(  m  f  (4  -f*)'£ 

For  a   small  oscillation,     ■typL>0  tityji.    —  Awl  a)t 
which  is  purely  sinusoidal.      As      ^        increases   it   is  distorted  as 
shown   in   Fig.    11(a). 

To  show  the  distortion  the  measure  A  is  found, 

A  -GWAf'+iSr)*  (40) 
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and  is  plotted  in  Fig„  10  along  f     „ 

(c)  Wave  form  of  f(t)   . 
From  Eq.  (37)  and  the  relation  (5), 

FffJ-fr^SfH       and       F'f^Sm 

rearranging   these   into  a  dimensionless   forms    we  get 

For  a   small  oscillation,     ^ „~Lk^    s    ^^ 

or"**  r" 

which  is  sinusoidal.  As  ^  increases  it  is  distorted  as  shown  in  Fig. 
11(b).  To  compare  this  wave  form  with  a  sinusoidal  one8  the  mea- 
sure B  is  found s 

which  is  plotted  along  ^  in  Fig.  10. 

(d)  Discussion. 

This  spring  is  of  softening  type,,  and  at|s  0„82  it  reaches 
its  maximum.   After  this  point  it  drops  downs  and  at  ft^fe      the  spring 
force  becomes  zero. 

The  measure  B  does  not  change  appreciably  up  to  j#s  Q„829 
and  Fig.  11(b)  shows  that  the  wave  form  is  not  distorted  much. 

For  any  value  of  p   which  indicates  the  non linearity  of  t 
spring,  the  above  solution  can  be  applied. 
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Eo   Example  5.   Piecewise- linear  softening  spring* 


Suppose  a  symmetrical  spring 
(as  shown)  is  described  by  the  equations, 

and  it  is  considered  the  case  of  ^S^only, 

The  equation  of  free  vibration  is 

n't  +  Sfr)  =*o 


(43) 


When  the  amplitude  of  deflection  ft    is  less  than  C  s  it  is  linear 9 
and  solution  will  not  be  included  here.   The  case  of  $}C   will  be  consider- 


ed. 

Let  us  denote  %,(t)  for 

fc(t)        for 
where  tc   -  the  time  when     yfc)  s  C 
and  let  uS   %f   and    lt£  -  ^6ff 

Then  Eq.  (43)  becomes 

i  +  ** =  * 

The  general  solutions  of  Eqs.  (43a)  are 


o  4  *  ^  tc 


o  4  t  £%. 


o 


^t^tc 


tw-n^^+ij+o-C)      £.±t^% 


(43a) 


(44) 


*  In  section  E  and  F  (which  follows),  for  convenience  in  notations 
only  the  first  quarter  cycle  of  motion  is  considered. 
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where  Q,    £)a,  <pt    ,   and  (p    ,  are  arbitrary  constants ,  which  together  with 
7c   may  be  determined  by  the  following  boundary  conditions; 

ty(c)  —  D  as  an  initial  condition. 

kr(p)   a*  -2  M(&)  by  the  conservation  of  energy. 

i<tj  -fttt) 

Solving  for  the  constants  in  Eqs.  (44),  we  get 

R  =  -«ci--jjf;  +  |  (45) 

To  obtain  the  frequency  response  curves,  the  period  of  free  oscil- 
lation (  P  =  0)   %  can  be  found: 

because  /fc(t)     reaches  its  maximum  value  of  #•  at  T -  %  „   Solving  for  J0 
we  get:    T.  =  -|  £J  -  &  ) 

^-flfcX'-*-) 

where   "I^*  -rr    is  the  period  of  a  small  free  oscillation. 

By  making  use  of  Eq.  (16)  for  forced  vibrations,  the  period  of  forced 
vibration  is  obtained, 

As  a  specific  example  let  us  take  tZt  —  A- -Ikx    in  the  above  symmetrical 
spring.   Then  we  get  (Jt  —2U)%     ,    and  the  solutions  of  Eqs.  (44)  become 

^         (44a) 
k(t)  =(f+3c)*«(*Ut+b)  -3c  tc^t6% 


where 


Xj=  [Ca  +  2£(j-c)+i^-c)aj 


£ 
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(a)  Frequency  response. 
From  Eq.  (46), 

where  /°  is  an  assignable  constant  as  defined  in  the  relation  (5).  In 
Fig.  12  the  frequency  response  curves  are  shown  for  f"  =  o  ,  ±o.sA,C  9 
and  dt  Ac  . 

(b)  Wave  form  of  wt)    • 

Eq.  (44a)  is  rewritten. 

For  g-^-C-     the  motion  is  purely  sinusoidal9  and  as  <*  in- 
creases it  is  distorted.  At  around  ^=-Z*5~c  i    distortion  reaches  a  maxi- 
mum;, and  for  further  increase  in  #•   the  wave  form  approaches  a  sinu- 
soidal form  again.   In  Fig.  13(a)  these  two  wave  forms  are  shown:   one  is 
sinusoidal  and  the  other  shows  the  maximum  distortion  at  ^  —2-$C  . 

To  show  the  distortion  the  measure  A  is  found9 

which  is  plotted  in  Fig.  12  along  %/q       .   Actually  the  distortion  is  s© 
small  that  the  measure  A  is  not  much  greater  than  unity8  where  the  maxi- 
mum measure  A  ■  1.023  at   ft—i-SC-     . 

(c)  Wave  form  of  f  (t) 

From  Eq.    (45)   and  the  relation    (5)j,    and  with  the   spring 
force      SW^iX+Mfa    ,\^C      >      for  this  spring,   we  get 
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*m.  04t4t 


F  J+3 

I*L  =    AMii&t+ba)  ic^r^%  (48) 

In  Fig.  13(b)  the  wave  form  is  plotted  for  f£C   ,  J  =Xc 
and  £c  .   Since  the  spring  force  has  a  slope  discontinuity  at  £  =  c  » 
the  wave  form  of  f(f)  is  slope  discontinuous  at  "t—t^ 

For  ^^  C   it  is  a  pure  sinusoidal  form.   As   ^   increases 
it  becomes  distorted;  maximum  distortion  occurs  at   $f=2>.$"c  approximately. 
Further  increases  in  H.     bring  out  a  sinusoidal  form  again*  similar  to 
$(~t)      wave  forms.   Every  wave  form  of  f(t)    does  not  go  beyond  a  limiting 
curve  which  is  marked  as  "Envelope"  in  Fig.  13(b). 

To  compare  the  wave  form  of  f(t)    with  a  sinusoidal  form, 
the  measure  B  is  found, 

8  -  -(SlrO-  <""&■)  jfa  +  W%-«*( « ffU  (W) 

which  is  plotted  along  /f/.  in  Fig.  12.   At  if--2-5c     8  approximate lya  3 
reaches  its  maximum  value  of  1.22.  At  this  point  the  distortion  is  a 
maximum. 

(d)  Discussion. 

The  solution  of  free  vibration  is  found  by  a  classical 
method,  and  for  the  forced  vibration  Harvey's  method  is  applied  to  that 
of  free  vibration. 

Fig.  13(b)  shows  that  J-(t)    is  not  far  away  from  a  pure 
sinusoidal  form*  and  the  measure  B  is  not  much  greater  than  unity. 
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F.   Example  6.   Piecewise  linear  hardening  spring 

Let  us  consider  a  symmetrical 
spring  as  shown   Let  4-K(  =  ft*   then  24/,«#4 

For  free  vibrations, 

Let  the  initial  conditions  be 

fie)  ~0     ,        *tf(0)  =2U(f) 
Then  the  solution  is  obtained  from  Eqs.  (44)  and  (45) , 

77         1 


where 


(a)   Frequency  response. 

From  Eq.  (46)  and  with  2U,\=u.L       ,  we  get 
X  -  .   (±-     &_)  '     - 


(51) 


In  Fig.  14  the  frequency  response  curves  are  plotted  for 
p   ■  Oj,   il,!!  ,  and  ±S$iC°      For  ^C  ,  it  is  linear  and  the  ratio 
(  T^Tn  )  is  unity.   And  as  H-   becomes  much  greater  than  C  ,  the  spring 

characteristic  again  approaches  linearity  and  the  ratio  (  ^e/'Tn      )  approaci 
0.5. 

(b)  Wave  form  of  ffi)   . 

Eqs.    (50)   are   rewritten: 

0  4tt  tc 


(52; 


In  Fig.  15(a)  the  wave  form  of  j}(t)     is  shown  for  ^^^ 


and  ^=/5'd'.   For  )f^C   »  the  spring  is  linear  and  the  wave  form 


purely 
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sinusoidal.  When  <l  becomes  greater  than  C     s  the  wave  form  is  distorted, 
at  approximately  %   ■  \-^C      the  distortion  reaches  a  maximum.   Fox  t 
increase  in  5   ,  the  spring  characteristic  approaches  linearity  again^, 
the  wave  form  approaches  a  sinusoidal  form. 

To  show  the  distortion  the  measure  A  is  found, 

which  is  plotted  along   fc£      in  Fig.    14.     At    ft^lsc  the  measure  A  reaches 
a  minimum   (A  =  0.965) . 

(c)     Wave   form  of  f[f)  . 

The   spring   characteristic   is 

5{;j)=Lf-lkc     for    yt\;?c 

Substituting  Eqs.  (50)  into  the  above  equation,  we  get 

f        n  ~iC 

Ff)    -    /unitLt+b)  tc^t-^%.  (54) 

In  Fig.  15(b)  the  wave  form  of  p(t)  is  shown  for  !f^C   8 
,  and  3.0  C  .   Tiere  is  a  limits  beyond  which  the  wave  form 
does  not  g   as  indicated  as  "Envelope"  in  th*  '-'jure. 

To  compare  the  f-(t)  wave  form  with  a  sinusoidal  form  the 
measure  B  is  founds 

B-sfep—fe)^*  +^W&Kk)  (55) 

which  is  plotted  in  Fig.  14  along  ££  .   For  ^<.c  the  spring  is  linear  and 
p(£j   is  a  pure  sinusoidal  form,  so  the  measure  B  is  unity.   As  H       goes 
beyond   c   ,  the  measure  B  decreases  rapidly >  and  at  ^ =I5C »  approximate lya 
B  a  0,7^5  as  a  minimum.   For  further  increase  in  H.        s    the  measure  B 
slowly  approaches  unity  where  the  spring  characteristic  becomes  linear  ag 
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(d)  Discussion,, 

This  example  is  similar  to  the  previous  one^  from  which 
the  general  solution  is  obtained.  The  measure  B  and  the  wave  form  of 
F  (t)     shown  in  Fig.  15(b)  tells  us  that  the  external  force  ("natural") 
is  not  80  far  away  from  a  sinusoidal  form. 


38 


40 


BIBLIOGRAPHY 

1.  T.  J,  Harvey ,  "Natural  Forcing  Functions  in  Nonlinear  Systems" , 
J.  Appl.  Mech.,  Vol.  25,  Trans.  ASME,  Vol.  80s  1958,  p  352. 

2.  M.  A.  Sadowsky8  "Nonlinear  Springs",  J.  Franklin  Inst.  Vol.  240 , 
No.  6,  Dec.  1945,  p  469. 

3.  U.  Ojalvo  and  G.  L.  Bleckman,  "On  Duffingss  Equation"9  J.  Appl. 
Mech.,  March  1961. 


41 


